Landau diamagnetism and magnetization of interacting diffusive conductors 
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We show how the orbital magnetization of an interacting disordered diffusive electron gas can be 
simply related to the magnetization of the non-interacting system having the same geometry. This 
result is applied to the persistent current of a mesoscopic ring and to the relation between Landau 
diamagnetism and the interaction correction to the magnetization of diffusive systems. The field 
dependence of this interaction contribution can be deduced directly from the de Haas-van Alphen 
oscillations of the free electron gas. Known results for the free orbital magnetism of finite systems 
can be used to derive the interaction contribution in the diffusive regime in various geometries. 



PACS Numbers: 

In recent years, there has been many theoretical works 
on the thermodynamic properties of mesoscopic elec- 
tronic systems, in particular concerning their orbital 
magnetism The simplest description of metals 

deals with non-interacting electrons in the absence of 
disorder. The correction to Landau susceptibility due 
to electron-electron interactions and phase coherence has 
been worked out by Altshuler et al. |2) . Similarly, the per- 
sistent current in mesoscopic rings has been extensively 
studied. The simplest description of this effect was first 
done in a strictly one-dimensional (ID) picture of non- 
interacting electrons and the effect of diffusion and 
interaction was described later by Ambegaokar and Eck- 
crn || and Schmid 

The very simple approach for free electrons and the 
more sophisticated description of interacting electrons in 
a disordered potential have been developed in a com- 
pletely independent way. Here, we show how these de- 
scriptions are closely related. The main result of this 
letter is a simple relation between the response of a clean 
non-interacting electron gas and the response of a diffu- 
sive electron system in the presence of interactions. This 
result originates from the very similar structures of the 
Schrodinger equation and of the diffusion equation which 
describe the two systems. 

As a first example, we show how the persistent cur- 
rent of a ID ballistic ring is related to the current of 
a quasi-lD diffusive ring in the presence of interactions 
0. Then we show how the interaction contribution to 
the orbital magnetism of any diffusive system can be de- 
duced immediately from the orbital response of the same 
non-interacting system. As a second example, we show 
how the interaction contribution to the susceptibility of a 
bulk diffusive system is derived directly from the Landau 
susceptibility. Then, from the de Haas-van Alphen oscil- 
lations of the free electron gas, we deduce the field depen- 
dence of the interaction induced magnetization. Finally 
we use this mapping to derive the finite size corrections 
(in L v /L) in the diffusive case from the 1/kpL correc- 
tions of the magnetization of the clean system. 

Classically, the probability p(r, r',u>) for a particle to 
diffuse from a point r to another point r' is solution of 



the diffusion equation 

(-^ + 7-L>V^)p 7 (r,r',w) = 5(r - r' 



(1) 



D is the diffusion coefficient. This probability has ac- 
tually two parts, a purely classical one (the Diffuson) 
and an interference part (the Cooperon) which results 
from interference between time reversed trajectories. The 
Cooperon has to been taken at r = r'. In a magnetic 
field, it obeys eq. (Q) where V has to be replaced by 

V + 2ieA/hc, A being the vector potential. The charge 
(— 2e) accounts for the pairing of time reversed trajecto- 
ries which are supposed to propagate coherently up to a 
time T^. 7 = and L v — U Dt v is the phase coher- 
ence length. 

The probability p(r, r',oj) has the same structure 
as the disordered averaged (retarded) Green's function 



G e (r,r',e) of the Schrodinger equation for a free parti- 
cle of energy e and charge — e in a disordered potential: 

( e -^ + ^)Sf(ry) = tf (,-0 (2) 

where r e is the elastic mean free path. In a field, V — > 

V + ieA/hc. The solutions of eqs.([l],||) can be written as 



Pj{ r , r ') = Pli v ; v '; 



0) 



^C(r)TMr') 



and 



G R (r',r,e) 



E 



(3) 



(4) 



where the eigenvalues E^' s are the solutions of similar 
equations 

*2 



DAib n = E*fa 



2m 



Aip n = E s J) n 



(5) 



with the mapping from the diffusion to the Schrodinger 
problem: 



D 

2e 
Try 



h 

2m 



2t> 



(6) 
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It has long been recognized that a Diffuson (or a 
Cooperon) behaves like a free particle with an effective 
mass m* = h/2D. The goal of this letter is to study the 
consequences of this mapping on the orbital magnetism 
of clean and diffusive systems. 

For a disordered finite system of size L, the Thouless 
energy E c , given by S.D/L 2 , is equivalent to the mean 
inter level spacing A = h 2 /2mL 2 of the eigenvalues of the 
Schrodinger equation. More interesting is the relation 
deduced from eq. (0). 



L 
I ■ 



-ikpL 



L 
2L 



(7) 



where l e = vpT e . l/r e spreads the levels of the 
Schrodinger equation while 1/t v spreads those of the 
diffusion equation. Inelastic disorder on the Cooperon 
plays thus the same role as elastic disorder on a free par- 
ticle. More important, the relation (^) expresses that 
the limit kpL ^S> 1 for the clean system corresponds to 
the macroscopic limit L 3> L v . Inversely the mesoscopic 
limit L <C L v corresponds to having only one Schrodinger 
particle in a box (kpL <C 1). 

Let us now apply this mapping to the calculation of 
the magnetization. First, the T = OK magnetic moment 
of the free electron gas (including spin) can be written as 



(8) 



where Af(e, B) is the double integral of the total density 
of states p(e, B). This contribution is known as the Lan- 
dau magnetization. Then taking into account electron- 
electron interactions in the Hartree-Fock picture gives an 
additional contribution . For a completely screened in- 
teraction U(r — r') = US(r — r') ||, this contribution is 
given by 



(M ei 



■755 / '"'M** 



-U— I (p{r,uj 1 )p(r,uj2)}drduj 1 duj2 (9) 



This expression contains the Hartree and Fock contribu- 
tions. n(r) is the local density. p(r, u>) is the local den- 
sity of states (per spin direction). The average product 
(p(r,ui)p(r,u>2)) is nothing but the Fourier transform of 
the return probability p y (r,r,t), so that one gets finally 
0: 



,„,.,, X h d f PJt) , 



(10) 



where P 7 (t) = / p 7 (r, r, t)dr in the space integrated re- 
turn probability. Ao = Upo is a dimensionless interaction 
parameter and po is the average density of states (per 
spin direction). Writing the density of states as 



p(e) = — — [ lraG R e(r, r)dr 

7T J 

and the integrated return probability as 



Pj(t)dt = / p 1 {v,v)dv 



(11) 



(12) 



one obtains immediately from eqs.(||,^0|) that the two 
magnetizations are related (since the 1/Pterm in cq.(|lO|) 
is equivalent to a double integral over 7) 

M=-^-Im[(M ee )( 7 =-^-iO)] (13) 

The sign = means that the two quantities are equal, pro- 
vided the substitutions (g) have been made. It should 
then be remembered that eq.(|lC|) corresponds to taking 
the first order contribution in Ao to the grand potential. 
It is known that taking into account higher diagrams in 
the Cooper channel, one has to renormalize the inter- 
action parameter which becomes energy dependent A(e) 

Ao^A(e)=A /(l + A m^) = l/m^ (14) 

where To is defined as To = epe 1 ^ . Then the relation 
( |l3| ) can be simply modified as: 



M = -lim Ao ^OT-Im[(M ee )( 7 
Ao 



iO)] (15) 



As an example, we consider the case of a ID diffusive 
ring of perimeter L pierced by a Aharonov-Bohm flux 
<f>. Starting from the flux dependent part of the return 
probability 



P(t) 



4nDt ^ 

p=i 



4Dt cosAnpip 



(16) 



where ip — (f)/4>o, (f>o being the flux quantum, one sim- 
ply gets from eq. (|Io|) , the harmonic dependence of the 
average persistent current due to interactions. 



16A ^V^(1 



V 



sin47rp< / 9 (17) 



This result, for L v = 00, was first obtained by Ambe- 
gaokar and Eckern (AE) ||. It was then generalized to 
the case where L v is finite fl2]] . Using the relation jl5|), 
one deduces immediately the average persistent current 
for a clean ID ring (clean means here that there is no 
diffusion. Disorder is only taken into account by a finite 
mean free path l e — vpr e ): 



I = l lQ y I( C0 Spfc F L - Sin P kFL ) e -pL/2le sin2w 

7T ^— ' p pkpL 
p—1 

(18) 



2 



with Iq — evp/L. This result has first been obtained for 
the case kpL 1 (in this case, the s'mx/x term cancels) 
in the absence of disorder (l e = oo) 0. Note that the 
correspondance between the AE current and the current 
of the ballistic ring is not trivial. The leading term in 
UfL for the clean case originates from the leading term 
in L/L^ in the diffusive case. Therefore, taking simply 
the AE result for the mesoscopic limit (L v — oo) would 
not have produced the correct result for the clean ring. 
In other words, the kpL 3> 1 limit corresponds to the 
macroscopic limit for the diffusive case. We will return 
to this point later where we will show how to derive the 
L v / L corrections to diffusive magnetization from perime- 
ter corrections in the ballistic case. 

Deducing the magnetization of a clean system from 
the one of the interacting system may not appear as the 
most useful procedure. More interesting is to deduce the 
properties of an interacting medium from those of the 
non-interacting one, i.e. to invert eq.(13). This inversion 
is given by: 



(M ee ) = 



M(e) 



de 



(19) 



with the substitution (J6j) . Defining M the magnetization 
of a free particule of mass h~/2D and charge 2e, so that 
M(e) = M(e), one can rewrite 



it Jo 



M(e) 
e + hj 



de 



(20) 



Again, recognizing that Cooper Channel renormalization 
modifies the interaction parameter, the energy depen- 
dence of this parameter can be incorporated exactly in 
the integral so that: 



1 



W e + fi7 



This is the main result of this paper. It gives straight- 
forwardly the magnetization of an interacting electron gas 
in terms of the magnetization of the same non-interacting 
system. 

As a example, we now consider the orbital response of 
a 2D clean system. The (spinless) Landau susceptibility 
gives the non oscillating part of the orbital response. It 
is given by x( e ) — ~ e 2 /(24irm) and is independent of 
the energy \{ e ) — XL- Then, using the mapping (||), 
the susceptibility of the cooperon is x(e) — ~if = 
— Axh{^FT e ) /h. From eq. J2l|) , one immediately deduces 
the interaction part of the susceptibility 



An ultraviolet cut-off l/r e has been added in order to 
cure the divergence at large energy. 

In 3D, the Landau susceptibility becomes energy de- 
pendent x(e) = — e 2 kp(e)/ (2Ai: 2 m) oc y/e, so is the sus- 
ceptibility x(e) of the cooperon, x(e) = — &XL-\f^- Con- 
trary to the 2D case where the susceptibility was constant 
in energy and of order ei?r e , integration in energy gives 
here a much smaller contribution. Using eq. (pi]) , one gets 
the interaction correction in 3D: 



Xe 



1G 



1 



\Xl\ ttV3 In T T e /h 



(23) 



Consider again the 2D clean case. In addition to the 
Landau contribution, the de Haas-van Alphen effect ex- 
presses the oscillatory behavior of the grand potential in 
1/B, with the fundamental period 1/Bq = ehjmtF- The 
grand potential is given by pjfl: 



SA(B) = -\xlB 2 



12 



2-kscf \ 
■ cos — J 



huj c 



(24) 



and the magnetic moment at fixed Fermi energy is given 
by M = —dSA/dB. Its dependence versus field has the 
well-known saw-toothed behavior. One may wonder how 
this behavior translates into the language of interacting 
diffusive electrons. To simplify, we restrict ourselves to 
the first order in Ao- Using the mapping ( |l9| , |20| ), one de- 
duces the interaction contribution to the magnetization, 
in units of XqUD/^: 



(M e . 



3 T R 



_8__5_ 
~^dB 



B 



E 



(- 



^/(2tt S — ) 



B 



(25) 



where f(x) = — c\(x) cos(x)— si(x) sin(a;). The fundamen- 
tal frequency Bq has been transformed into the charac- 
(21) teristic field B v = %/ (AeDr^). Alternatively, this mag- 
netization could have been obtained from eq. (|l0|) , with 
the following expansion of the return probability in a 
constant field: 



P(t) = 



sinh An BDt/<j) 
1 



4nDt 



1 



t z 



t 2 + a 2 s 2 



(26) 



with a — 4>q/ (ABD) . The result ( p5| ) can be easily gen- 
eralized to all orders in Ao by considering the explicit 
dependence A(e) is eq.(f2l"l). 



AhD , ]nT T v /h 
In — — 

3 4>l \nT T e /h 



h lllToT e /?i 



(22) 
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FIG. 1. Magnetization of a diffusive interacting electron 
gas calculated to first order in Ao, in units of \ohD/4>i- The 
dashed line shows the linear low field behavior, see eq. 
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FIG. 2. Susceptibility of a diffusive interacting electron 
gas in units of XohD/<j>Q. The amplitude at zero field is 
4/31nr v ,/r e , see eq. (E5h. 



Let us finally note that the limit kpL 3> 1 for 
the Schodinger equation corresponds to the macroscopic 
regime L L v for the diffusion equation. The oppo- 
site, so called mesoscopic regime L -C L v would corre- 
spond to kpL -C 1, for which only the ground state is 
occupied. In the diffusive context, this ground state is 
called the zero mode. The cross-over between the meso- 
scopic regime where only a few modes are relevant to the 
macroscopic regime where there is a quasi-continuum of 
diffusion modes is quite difficult to describe [^4[ . It is 
then quite useful to know the finite size 1/kpL correc- 
tions to the Landau susceptibility which have been ex- 
tensively studied |^5|. These corrections are usually of 
the form |15|: 



X (L) ~ X (oo) 1 



(27) 



Thus, knowing the finite size corrections to the Landau 
diamagnetism, one can get the L v / ' L corrections to the 
bulk susceptibility Xee ■ For L^> L v , they are of the form 



(28) 



In conclusion, we have shown that the magnetization of 
a diffusive interacting electron gas can be deduced from 
the magnetization of the non-interacting system. This 
mapping allows the study of finite size properties of dif- 
fusive systems, in particular the cross-over between the 
macroscopic and the mesoscopic regimes. 
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